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ABSTRACT 

The high-frequency-peaked BL Lac (HBL) lES 0806+524 (z = 0.138) was discovered in VHE y rays in 2008. 
Until now, the broad-band spectrum of lES 0806+524 has been only poorly characterized, in particular at high 
energies. We analysed multiwavelength observations from y rays to radio performed from 2011 January to March, 
which were triggered by the high activity detected at optical frequencies. These observations constitute the most 
precise determination of the broad-band emission of lES 0806+524 to date. The stereoscopic MAGIC observations 
yielded a y-ray signal above 250 GeV of (3.7 ± 0.7) per cent of the Crab Nebula flux with a statistical significance 
of 9.9 <T. The multiwavelength observations showed significant variability in essentially all energy bands, including 
a VHE y-ray flare that lasted less than one night, which provided unprecedented evidence for short-term variability 
in lES 0806+524. The spectrum of this flai'e is well described by a power law with a photon index of 2.97 ± 0.29 
between ~150 GeV and 1 TeV and an integral flux of (9.3 ± 1.9) per cent of the Crab Nebula flux above 250 GeV. The 
spectrum during the non-flaring VHE activity is compatible with the only available VHE observation performed in 
2008 with VERITAS when the source was in a low optical state. The broad-band spectral energy distribution can be 
described with a one-zone Synchrotron Self Compton model with parameters typical for HBLs, indicating that lES 
0806+524 is not substantially different from the HBLs previously detected. 

Key words: BL Lacertae objects; individual: lES 0806+524 - galaxies: active - galaxies; jets 
- gamma rays: galaxies - radiation mechanisms: non-thermal 


1 INTRODUCTION 

Active galactic nuclei (AGN) are among the most variable objects 
in the known Universe. Their broad-band energy spectrum span¬ 
ning from radio to very high energy (VHE, E > 100 GeV) y rays 
can be characterised by two distinct peaks: one in the sub-mm to 
X-ray range, commonly interpreted as synchrotron radiation, and 
a second one in the y-ray band that is hypothesised to originate 
from inverse Compton scattering of photons. The emission is as¬ 
sumed to originate from relativistic particle jets launched along 

* CoiTesponding authors: C. Schultz, E-mail: comelia.schultz@pd.infn.it, 
S. Buson, E-mail: buson@pd.infn.it, R. Reinthal, E-mail: rirein@utu.fi, F. 
Tavecchio, E-mail: fabrizio.tavecchio@brera.inaf.it 


the axis of the accretion disc of matter that surrounds a supermas- 
sive black hole (e.g. Blandford & Rees 1974). Numerous AGNs 
are known to be blazars, characterised by relativistic jets closely 
aligned to the line of sight of the observer. Based on their optical 
spectra, blazars are divided into two classes: flat spectrum radio 
quasars (FSRQ) that show broad emission lines, and BL Lacer¬ 
tae objects (BL Lacs) characterised by the weakness or even ab¬ 
sence of such emission lines (Weymann et al. 1991; Stickel et al. 
1991). Depending on the frequency of the low-energy peak of the 
spectral energy distribution (SED) the latter class is subdivided 
into high- (HBL), intermediate- (IBL), and low- (LBL) frequency- 
peaking BL Lac objects (Padovani & Giommi 1995). The Tuorla 
Blazar Monitoring Program (Takalo et al. 2008) monitors a sample 
of objects in the optical band, triggering VHE observations by the 
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MAGIC (Major Atmospheric Gamma-Ray Imaging Cherenkov) 
telescopes during high optical states (Albert et al. 2006, 2007a, 
2008a; Anderhub et al. 2009; Aleksic et al. 2012a,b). This allowed 
the discovery in VHE y rays of several blazars, especially HBLs 
and leads to the suggestion that in HBLs a high optical state is typ¬ 
ically accompanied by a high state in VHE y rays. In some cases, 
studies of broad-band emission confirm a possible connection be¬ 
tween the two wavebands (e.g. Aharonian et al. 2009; Aleksic et al. 
2012c) whereas in other sources this seems to be only partially the 
case (e.g. Foschini et al. 2007, 2008). 

lES 0806-H524 (RA 08:09:49.18673, DEC 
52:18:58.2507; J2000, Petrov & Taylor 2011) is classi¬ 

fied as a BE Lac (Schachter et al. 1993) with a redshift of 
0.138 (Badeetal. 1998). It was suggested as a VHE candidate 
(Costamante & Ghisellini 2002) with a predicted intrinsic flux 
of F£>o. 3 TeV = 1-36 X 10“*' cm“^ s“'. Several VHE observa¬ 
tions have been carried out by the Whipple Collaboration and 
the HEGRA Collaboration yielding flux upper limits above 
300 GeV (Horan et al. 2004; de la Calle Perez et al. 2003) and 
1.09 TeV (Aharonian et al. 2004), respectively. 

In 2008, the VERITAS Collaboration reported the first de¬ 
tection of lES 08064-524 in the VHE y-ray band (Acciari et al. 
2009). The collected data set spanned from 2006 November to 
2008 April yielding a total of 245 excess events with a signifi¬ 
cance of 6.3 (T. The integral flux above 300 GeV corresponded to 
1.8 per cent of the Crab Nebula flux, which is below earlier up¬ 
per limits obtained by MAGIC-I observations in 2005 (5.6 per 
cent C.U.' above 230 GeV; Albert etal. 2008b; 7.2 per cent C.U. 
above 140 GeV; Aleksic et al. 2011). No significant variability on 
a monthly timescale (the for a fit with a constant is 

6.78/5; Acciari et al. 2009) could be established for this object. 
The spectrum was obtained only for a subset of data taken with 
four telescopes during winter 2007/2008. Accordingly, the spec¬ 
tral characteristics were poorly defined. Thus, the spectral index 
was measured to be 3.6 ± L0stat±0.3sys in a narrow range between 
300 GeV and 700 GeV. The time-independent version of the one- 
zone jet radiation transfer code of Bottcher & Chiang (2002), with 
parameters appropriate for a pure Synchrotron Self Compton (SSC) 
model (Maraschi & Tavecchio 2003) was able to describe the data 
sufficiently well within the uncertainties. 

In this paper we report the highly significant detection of this 
source with the MAGIC stereoscopic system during a flaring state 
in 2011 February. We analyse in detail the variability and the spec¬ 
tral evolution of the source and present complementary data from 
observations in high-energy (HE, 30MeV < E < 100 GeV) y-rays 
carried out by the Fermi Large Area Telescope (LAT), in X-rays 
performed by the Swift satellite and in the optical R-band by the 
Kungliga Vetenskapsakademien (KVA) telescope. Radio data cov¬ 
erage is provided by the Owens Valley Radio Observatory (OVRO) 
telescope at 15 GHz. Thanks to the good multi wavelength (MWL) 
coverage, we investigate the connection among these wavebands. 
Finally, we model the SED including all MWL data. 


' The integrated flux level in Crab units (C.U.) is obtained by normalizing 
the integrated flux measured above a certain threshold to the Crab Nebula 
flux, which is considered to be stable, measured above the same threshold. 
^ d.o.f.= Degrees of freedom. 


2 OBSERVATION AND DATA ANALYSIS 

2.1 MAGIC observations and data analysis 

Since 2009 MAGIC has operated as a stereoscopic system of two 
17 m Imaging Atmospheric Cherenkov Telescopes (lACTs) located 
at the Roque de Los Muchachos, La Palma, Canary Islands (28.8° 
N, 17.8° W, 2225 m a.s.l.). Due to its low energy threshold (as low 
as 60 GeV in normal trigger mode) and high sensitivity^, it is a well 
suited instrument for VHE y-ray observations of blazars. 

MAGIC observed lES 08064-524 between 2011 January and 
March on 13 nights for a total exposure of ~24h. The observa¬ 
tions from Eebruary to March were triggered by a flux increase in 
the optical R-band as part of a dedicated MAGIC Target of Oppor¬ 
tunity (ToO) program to observe blazars showing high activity at 
other wavelengths. After applying quality selection cuts based on 
the event rate, ~3.8h of data were discarded. Additionally, correc¬ 
tions for the dead time of the readout system yielded an effective 
observation time of 16.1 h. Part of the observations were carried out 
under moderate moonlight conditions. 

Observations were performed in the so-called wobble 
mode (Fomin et al. 1994) during which both telescopes alternated 
every 20 minutes between two sky positions with an offset of 0.4° 
from the source. The acquired data cover a zenith angle from 23° 
to 46°. 

The data analysis was performed using the MAGIC stan¬ 
dard tool ‘MARS’ (Morale)o et al. 2009) including adaptations 
to stereoscopic observations. Based on the timing informa¬ 
tion (Aliu et al. 2009), and absolute cleaning levels, the image 
cleaning was performed. The images were parametrized in each 
telescope individually according to the description of Hillas (Hillas 
1985). 

For the reconstruction of the shower arrival direction the ran¬ 
dom forest regression method (RF DISP method; Aleksic et al. 
2010) with the implementation of stereoscopic parameters such as 
the impact distance of the shower on the ground and the height of 
the shower maximum was used (Lombardi et al. 2011). 

To perform the gamma-hadron separation, the random forest 
method was applied (Albert et al. 2008c) using the image parame¬ 
ters of both telescopes and the shower impact point and height max¬ 
imum. Energy lookup tables were used for the energy reconstruc¬ 
tion. An angular resolution of ~0.07° at 300 GeV and an energy 
resolution as good as 16 per cent in the medium energy range (few 
hundred GeV) are achieved (details on the stereo MAGIC analysis 
can be found in Aleksic et al. 2012d). 

For sources with VHE y-ray spectra similar to that of the Crab 
Nebula, the sensitivity of the MAGIC stereo system is best above 
250 - 300 GeV. Eor sources with spectral shapes softer than that 
of the Crab Nebula, the best performance occurs at slightly lower 
energies. Consequently, we chose 250 GeV as the minimum energy 
to report signal significances and y-ray fluxes in light curves, while 
for the spectral analysis, in order to use all the available informa¬ 
tion, we also considered energies well below 250 GeV, where the 
analysis of the MAGIC data can still be performed. 

2.2 Fermi-LAT data analysis 

lES 08064-524 has been observed by the pair conversion telescope 
Fermi-LAT optimised for energies from 20 MeV up to energies be- 

^ Better than 0.8 per cent of the Crab Nebula flux in 50 h of observing time 
above 290 GeV (Aleksic et al. 2012d). 
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yond 300 GeV (Atwood et al. 2009). In survey mode the Fermi- 
LAT scans the entire sky every three hours. The data sample, which 
consists of observations between 2010 November 22 and 2011 June 
13, was analysed with the standard analysis tool gtlike, part of the 
Fermi Science Tools software package (version 09-27-01) available 
from the Fermi Science Support Center (FSSC)''. 

We selected P7CLEAN events located in a circular region of 
interest of 10° radius centred on the position of lES 0806+524. 
To reduce the contamination from the Earth-limb y rays produced 
by cosmic-ray interactions with the upper atmosphere, data were 
restricted to a maximum zenith angle of 100°. 

Eor the y-ray signal extraction, the background model used 
included two components: Galactic diffuse and an isotropic diffuse 
emission, provided by the publicly available files gaL2yearp7v6_- 
trim_v0.fits and iso_p7v6clean.txt^. 

The model of the region of interest (ROI) also included 
sources from the second Fermi-LPCt catalog (2EGL; Nolan et al. 
2012) that are located within 15° of lES 0806+524. These sources, 
as well as the source of interest, were modelled with a power law 
spectral shape with the initial parameters set to their 2FGL val¬ 
ues. In the 2FGL catalog lES 0806+524 is associated with the 
source 2FGL J0809.8+5218, which has been reported with a flux 
of (2.2 ± 0.3) X 10“® cm“^ s“* and a photon index of (1.94 ± 0.06) in 
the 2FGL catalog. When fitting, the spectral parameters of sources 
within 10° from our target were allowed to vary while those within 
10° - 15° were fixed to their initial values. 

During the spectral fitting, the normalizations of the back¬ 
ground models were allowed to vary freely. Spectral parameters 
were estimated from 300 MeV to 300 GeV using an unbinned max¬ 
imum likelihood technique (Mattox et al. 1996) taking into account 
the post-launch instrument response functions (IRF; specifically 
P7CLEAN_V6, Ackermann et al. 2012). 

During the MAGIC observing period, the source was not sig¬ 
nificantly detected on a daily basis. To ensure a good compromise 
between having a significant detection in most of the intervals and 
details on the temporal behavior of the source, the light curve was 
produced with weekly binning. Flux upper limits at 95 per cent 
confidence level were calculated for each time bin where the test 
statistic (TS®) value for the source was below 9 (see section 3.2). 

The systematic uncertainty in the flux is dominated by the sys¬ 
tematic uncertainty in the effective area, which is estimated to be 
10% at 100 MeV, decreasing to 5% at 560 MeV, and increasing to 
10% at 10 GeV (Ackermann et al. 2012). The systematic uncertain¬ 
ties are smaller than the statistical uncertainties of the data points 
in the light curve and spectra. 


2.3 Swift observations and analysis 

Beside its prime objective of detecting and following up y-ray 
bursts, since its launch in 2004 November the Swift Gamma- 
Ray Burst observatory has become an instrument suitable for 


^ http://fermi.gsfc.nasa.gov/ssc/ 

^ http://fermi.gsfc.nasa.gov/ssc/data/access/lat/BackgroundModels.html 
® The test statistic value quantifies the probability of having a pointlike y- 
ray source at the location specified. It corresponds roughly to the standard 
deviation squared assuming one degree of freedom (Mattox et al. 1996). 
The TS is defined as —21og(Lo/L), where Lq is the maximum likelihood 
value for a model without an additional source (i.e. the ’null hypothesis’) 
and L is the maximum likelihood value for a model with the additional 
source at the specified location. 


various purposes due to its fast response and its MWL abil¬ 
ities (Gehrels et al. 2004). Swift hosts three telescopes opti¬ 
mised for different energy ranges: the Burst Alert Telescope 
(BAT; Barthelmy et al. 2005) suited for observations between 
15-150 keV, the X-ray telescope (XRT; Burrows et al. 2005) 
with a 0.3-10keV coverage and the UV/Optical Telescope 
(UVOT; Roming et al. 2005) optimal for observations within the 
1800-6000 A wavelength range. 

Following the VHE y-ray flare detection of lES 0806+524 
by MAGIC (Mariotti et al. 2011), Swift ToO observations were re¬ 
quested and performed from 2011 February 26 to March 2; five ob¬ 
servations together with MAGIC. A high-activity state of the source 
was confirmed reporting clear variability in X-rays (Stamerra et al. 
2011). The source was monitored with the 5w//i/XRT in photon 
counting (PC) mode with ~2ks snapshots each night for a total ex¬ 
posure time of lOks. 

The data processing was performed with the xrtpipeline 
vO.12.6 distributed by HEASARC as part of the HEASoft pack¬ 
age. Events with grades 0-12 (according to the Swift nomencla¬ 
ture; Burrows et al. 2005) were selected for the PC data and the 
response matrices included in the Swift CALDB’ were applied. All 
observations showed a source count rate > 0.5 counts s“*, thus 
pile-up correction was required. We extracted the source events 
from an annular region with an inner radius of 3 pixels (estimated 
by means of the point spread function (PSF) fitting technique, 
see Moretti et al. 2005) and an outer radius of 30 pixels, while 
background events were extracted from an annular region centred 
on the source with radii of 70 and 120 pixels. Ancillary response 
files were generated with xrtmkarf, and account for different ex¬ 
traction regions, vignetting and PSP corrections. 

We fit the spectra in the 0.3-10 keV energy range. The spec¬ 
tra were rebinned with a minimum of 20 counts per energy 
bin so that the )d minimization fitting technique could be used. 
The spectral analysis was performed with XSPEC (Amaud 1996) 
adopting a simple power law model that includes a hydrogen- 
equivalent column density fixed to the Galactic value nFl = 4.1 x 
10“° cm-2 (Kalberla et al. 2005). 

Swi/r/UVOT data were taken with the ‘filter of the day’ (ei¬ 
ther U, UVWl, UVM2, or UVW2 filter; Poole et al. 2008) chosen 
day by day by the Swift science planners. We analysed the data 
using the uvotsource task included in the HEASOFT package. 
Source counts were extracted from a circular region of 5 arcsec 
radius centred on the source, while background counts were de¬ 
rived from a circular region of 10 arcsec radius in the source neigh¬ 
borhood. Conversion of magnitudes into dereddened flux densities 
was obtained by adopting the extinction value E{B — V) = 0.039 
from Schlegel, Finkbeiner & Davis (1998), the mean Galactic ex¬ 
tinction curve in Fitzpatrick (1999) and the magnitude-flux calibra¬ 
tions by Bessell, Castelli & Plez (1998). 

2.4 Optical observations and data analysis 

The optical data used in this study were obtained with the KVA* 
located at the Roque de los Muchachos observatory on La Palma, 
and consisting of two telescopes mounted on the same fork. The 
smaller, a 35 cm Celestron, is used for photometric measurements, 
while the larger, a 60 cm one, is used for polarimetric observa¬ 
tions. They are operated remotely from Finland. The photometric 

’ http://heasarc.gsfc.nasa.gov/docs/heasarc/caldb/swi£t/ 

* http: //WWW . astro. utu. fi/telescopes/6Slapalma. htm 
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measurements are performed in the optical R-band using differ¬ 
ential photometry, by having the target and the calibrated com¬ 
parison stars in the same CCD images (Fiorucci, Tosti & Rizzi 
1998). The magnitudes of the source and comparison stars are mea¬ 
sured using aperture photometry and converted to linear flux den¬ 
sities according to the formula FfJy) = Fg x where 

Fo is a filter-dependent zero point (Fo = 3080 Jy in the R-band, 
from Bessell 1979). To obtain the AGN core emission, contribu¬ 
tions from the host galaxy and possible nearby stars that add to the 
overall flux have to be subtracted from the obtained value. These 
values were determined by Nilsson et al. (2007) and in the case of 
lES 0806+524 the host galaxy contribution of (0.69 ± 0.04) mjy 
was subtracted from the overall flux. 

The KVA telescope is operated under the Tuorla Blazar Mon¬ 
itoring Program®, which has been running as a support program 
to the MAGIC observations since the end of 2002. The project 
uses the Tuorla 1 m (located in Finland) and the KVA telescopes 
to monitor candidates (from Costamante & Ghisellini 2002) and 
known TeV blazars in the optical waveband and to provide alerts to 
MAGIC on high states of these objects in order to trigger follow¬ 
up VHE observations. lES 0806+524 was one of the objects on the 
original target list and has therefore been monitored regularly since 
the beginning of the program. The object had been relatively dor¬ 
mant over the last years, showing some variability but no particular 
flaring activity. The large flare occurring at the end of 2010 (see 
section 3.2) is by far the largest optical activity recorded from this 
source during the last ten years, although this kind of flaring ac¬ 
tivity is not uncommon in the optical light curves of other blazars 
observed within the Tuorla monitoring program. 

2.5 OVRO data analysis 

Observations at 15 GHz were carried out with the 40 m OVRO 
radio telescope, located in California. The observations of lES 
0806+524 were carried out in the framework of a blazar monitor¬ 
ing program (Richards et al. 2011) measuring the source flux den¬ 
sity twice a week. Occasional gaps in the data sample are due to 
poor weather conditions or maintenance. Observations were per¬ 
formed using a Dicke-switched dual-beam system, with a second 
level of switching in azimuth to alternate between source and sky 
in each of the two horns, which removes much of the atmospheric 
and ground interference (Readhead et al. 1989). The data were cal¬ 
ibrated against 3C 286 with an assumed flux density of 3.44 Jy 
at 15 GHz (Baars et al. 1977). The data were analysed using the 
pipeline described in Richards et al. (2011). 

3 RESULTS 
3.1 MAGIC results 

Above an energy threshold of 250 GeV the MAGIC data yield (after 
event selection cuts) an excess of 148 ± 17 events in the distribution 
of the squared angular distance 9^ between the reconstructed event 
direction and the catalog position of lES 0806+524, i.e. in the so- 
called ‘On’ region. The background level of 114 ± 6 events was 
calculated from three equivalent ‘Off’ regions, located at 90“’, 180“’ 
and 270“’ with respect to the reconstructed source position in the 
camera, applying the same event selection cuts (Eig. 1). 

® Project web page: http://users.utu.fi/kani/lm/ 



Figure 1. Distribution of the squared angular distance (0^) for the on-source 
counts in the direction of lES 0806+524 (black points with error bars) and 
the normalised off-source events (gray histogram and open black squai'es) 
extracted from three background regions which are located at 90“’, 180“’ and 
270° with respect to the reconstructed source position in the camera. The 
signal is extracted in the 6^-region indicated by the vertical dashed line. 



Figure 2. VHE y-ray flux (E > 250 GeV) of lES 0806+524 for the single 
MAGIC observations performed in 2011. The red arrows correspond to 95 
per cent confidence upper limits, which were computed for the observations 
where the interval flux ± error contains zero. The measured fluxes for these 
intervals ai'e also shown (gray circles). 

The significance of the event excess corresponds to 9.9 cr cal¬ 
culated with formula [17] of Li & Ma (1983). During the MAGIC 
observations, the source underwent a flaring event on February 
24 (Mariotti et al. 2011). Within 3.0 h of observation an excess of 
50 ± 8 events (15 ± 2 background events) above 250 GeV was mea¬ 
sured corresponding to a confidence level of 7.6 cr. After excluding 
the flare of February 24 from the data set, the remaining MAGIC 
observations (13.1 h) still show a significant detection of 96 ± 15 
excess events (99 ± 6 background events) above 250 GeV corre¬ 
sponding to 7.3 <T. In the following, we will refer to the VHE y-ray 
flare as the high state. The low state refers to a low activity in VHE 
y rays including the remaining MAGIC observations. 

The integral flux above 250GeV was (5.9 ± 1.1) x 
10“'^ cm“^ s“' corresponding to (3.7 ± 0.7) per cent C.U. Exclud¬ 
ing the flare, the integral flux was measured to be (3.1 ± 1.0) x 
10“'^ cm“^ s“' corresponding to (1.9 ± 0.7) per cent C.U. (2.1 ± 0.7 
per cent C.U. above 300 GeV). This result is in good agreement 
with the flux level above 300 GeV reported by VERITAS in 2009 
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Figure 3. Unfolded high-state differential energy spectrum of lES 
0806-t524 measured by MAGIC on 2011 February 24. The open black 
triangles correspond to the measured spechum to which a simple power 
law (solid black line) is fitted; the filled blue triangles depict the mea¬ 
sured spectrum after correction for the EBL attenuation using the model 
by Dominguez et al. (2011). For comparison, the observed spectrum and 
derived upper limits (gray squares and an'ows respectively) published by 
VERITAS (Acciari et al. 2009) and the Crab Nebula spectrum (red dashed 
line) published by MAGIC (Albert et al. 2008d) are shown. See text for 
further details. 


Figure 4. Unfolded low-state differential energy spectrum of lES 
0806-t524 measured by MAGIC in 2011. The open black triangles 
correspond to the measured spectrum to which a simple power law 
(solid black line) is fitted; the filled blue friangles depict the mea¬ 
sured spectrum after correction for the EBL attenuation using the model 
by Dominguez et al. (2011). For comparison, the observed spectrum and 
derived upper limits (gray squares and arrows respectively) published by 
VERITAS (Acciari et al. 2009) and the Crab Nebula spectrum (red dashed 
line) published by MAGIC (Albert et al. 2008d) ai'e shown. 


(1.8 per cent C.U.; Acciari et al. 2009). During the night when the 
source was flaring the integral flux was (1.5 ± 0.3) x 10“'* cm“^ s“' 
equal to (9.3 ± 1.9) per cent C.U. The flare showed a flux of ~3 cr 
above the mean flux level'" which corresponded to a flux increase 
of about a factor of three (Fig. 2). No intra-night variability was 
found during the VHE y-ray flare within the statistical and system¬ 
atic uncertainties. Besides the first flare observed by MAGIC on 
2011 February 24 (MJD = 55616), the VHE data showed a hint of 
increasing flux towards March 2 (MJD = 55622), when the source 
was detected at ~5 cr with a flux of (7.0 ± 2.3) per cent C.U. above 
250 GeV. Fitting the overall light curve with a constant function 
yields a probability of 0.5 per cent ix^lA.o.f. = 24.92/10) for a non¬ 
variable source. Since the nights before and after the flaring event 
of February 24 showed a significantly lower flux, we assume short¬ 
term variability of the timescale of one day as an upper limit. 

The differential energy spectra of the high and low states 
are shown in Figs. 3 and 4 respectively. Both spectra can be 
described with a simple power law" of the form dNjdE = 
/o(£/500GeV)“*^ with a flux normalization /o at 500 GeV of 
(1.25±0.21)xl0“" cm“- 8“* TeV“' for the high and (4.39± 1.04)x 
10“*^ cm“^ s“* TeV“* for the low state. The photon index F was 
found to be 2.97 ± 0.29 for the high and 2.65 ± 0.36 for the low 
state. The spectra have been unfolded using the Tikhonov algorithm 
to correct for the finite energy resolution of MAGIC. Different un¬ 
folding algorithms as described in Albert et al. (2007b) were com¬ 
pared and found to agree within errors. 

To account for Extragalactic Background Light (EBL) attenu- 


*" The deviation from the constant flux level was computed by dividing the 
difference of the flux measured during the VHE y-ray flare and the mean 
flux level by the sum in quadrature of the statistical errors on the flux in the 
high and mean flux levels. 

" The fits of the high and low state spectra with a power law have a 
X^ld.o.f. of 0.88/3 (83 %) and 1.16/3 (76 %), respectively. 


ation, a correction by means of the EBL model of Dominguez et al. 
(2011)*^ was applied. The corrected spectral flux of the two dif¬ 
ferent source states follows a power law parametrized by a pho¬ 
ton index F = 2.30 ± 0.52 and 2.15 ± 0.59 as well as a flux nor¬ 
malization /o at 500 GeV (2.57 ± 0.81) x 10“" cm“^ s“* TeV“* and 
(9.10 ± 4.03) X 10“'^ cm“^ s“' TeV“*, for the high (Fig. 3) and low 
state (Fig. 4), respectively. The low state spectrum observed by 
MAGIC shows a good agreement within errors with the observed 
spectral points published by VERITAS. A comparison of the differ¬ 
ential energy spectra describing the different source states indicates 
an increase in flux but no significant hardening of the spectrum. 


3.2 Multiwavelength results 

Figure 5 shows the long-term MWL light curves in VHE y rays 
(MAGIC), HE y rays (Fermz-LAT), X-rays (Swift), in the R-band 
(KVA, Tuorla) and radio regime (OVRO telescope). 

The Fermi-LAI light curve was produced fixing the power law 
index of the source of interest to the value derived from the spectral 
analysis of the data set from 2010 November 22 to 2011 June 13, 
i.e. 1.88 ±0.17. 

Possible variations in the source emission in HE y rays have 
been tested following the same likelihood method described in the 
2FGL (Nolan et al. 2012). The method, applied to the data from 
2010 November to 2011 June, indicates that the source is vari¬ 
able (TS™r = 52 for 29 d.o.f.). When applied only to data taken 
during the MAGIC observation period (MJD = 55603 to MJD = 
55622) the probability for constant emission increases to 23 per 
cent (TS^r/d-O-f- = 10.5/8). Considering the entire data sample, i.e. 


*^ This model is among the most recent ones 
(e.g. Franceschini, Rodighiero &Vaccari 2008; Finke, Razzaque &Dermer 
2010; Gilmore et al. 2012) and agrees within the uncertainties with other 
models. 

























6 J. Aleksic et. al. 



Figure 5. MWL light curves of !ES 0806+524 from 2010 December to 2011 June. From top to bottom: VHE y rays (red circles) by MAGIC, HE y rays 
(orange triangles) by Fermi-LAT, X-rays (green squares) by Swift/XRT, R-band (purple triangles) by KVA, and radio (blue diamonds) by OVRO. The optical 
data are corrected for the host galaxy contribution according to Nilsson et al. (2007). Upper limits at 95 per cent conhdence level are indicated by downward 
arrows. In the VHE light curve, upper limits (red arrows) are computed for the observations where the interval flux ± error contains zero. The measured fluxes 
for these intervals are also shown (gray circles). The individual light curves are daily binned except for the light curve in the HE band from Fermi, where a 
binning of seven days has been used. The horizontal dashed lines report the mean flux in each light curve. 


from 2010 November 22 to 2011 June 13, the average flux level 
corresponds to (2.10 ± 0.14) x 10“* cm“^ s“' (E > 300 MeV). 

A smooth flux increase from the beginning of 2011 March 
until 2011 mid-April reaching a maximum of (5.2 ± 1.3) x 
10“* cm“^ s“* was observed with a delay compared to the other 
wavelengths. Given the long integration time of seven days, no 
clear conclusion regarding simultaneous source variability with re¬ 
spect to that observed in VHE y rays can be drawn within the time 
interval of the flare. In particular, during the night of the VHE y- 
ray flare the source was only marginally detected (TS = 4) in the 
Ferm/-LAT energy band. 

In order to provide simultaneous coverage in this energy band, 
we computed the 95 per cent confidence level upper limit be¬ 
tween 300 MeV and 300 GeV for February 24 corresponding to 
7.2 X 10“'* ergcm“^ s“* at a center energy of 580 MeV. 

The X-ray spectra derived with the individual 5wi/f/XRT ob¬ 
servations are reported in Table 1 . The time-averaged X-ray flux in 
the band 2-10 keV is (5.1 ±0.2) x 10“*^ ergcm“^ s“*, with a flux in¬ 
crease by about a factor of six when comparing the lowest flux (oc¬ 


curring on February 1) with the highest flux (occurring on March 
2). A fit with a constant function yields ^^/d.o.f. = 150/5, showing 
a clear X-ray variability and confirming the high-activity state of 
the source during the MAGIC observations performed at the end of 
February and beginning of March. 

SwiftfUNOT observations were carried out with different fil¬ 
ters in the ultraviolet bands. The brightness of (14.4 ± 0.03) mag 
measured on March 2 in the UVW2 and UVM2 bands is al¬ 
most unchanged with respect to February 1 where a brightness of 
(14.5 ± 0.03) mag was measured. However, lES 0806+524 appears 
about 1 mag brighter compared to the UV flux observed in 2008 
March. During the observations, the UV band photometry is com¬ 
patible with a constant flux within the errors. 

In the R-band, the core flux showed a large increase over 
the long-term base level starting from 2010 November (Fig. 6). In 
the subsequent months the flux density continued to increase until 
reaching a maximum of (4.75 ± 0.09) mjy (host galaxy subtracted) 
in the night of 2011 February 24 during the outburst in VHE y 
rays; this was almost three times higher than the quiescent state 
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Figure 6. Host-galaxy-corrected (Nilsson et al. 2007) long-term light curve of lES 0806+524 of the optical R-band observed by the KVA telescope. The gray 
shaded area marks the period of the MAGIC observations in 2011. The quiescent state (Reinthal et al. 2012) is indicated as a dashed line. 


Table 1. Log and fitting results of SwiftfKKt observations of lES 0806+524 
using a power-law model with Ah fixed to Galactic absorption in the range 
from 0.3 to lOkeV. 


Observation 
Date [MJD] 

Photon index 

r 

Flux 2-10 keV 
xl0“*^ erg 

y^/d.o.f. 

55593 

2.75 ± 0.30 

1.76 + 0.43 

Cash" 

55617 

2.36 ± 0.09 

4.91+0.36 

36.84/34 

55619 

2.44 + 0.10 

8.13 + 0.73 

27.73/27 

55620 

2.38 ± 0.09 

5.73 + 0.51 

46.39/37 

55621 

2.40 + 0.11 

5.63 + 0.49 

33.67/27 

55622 

2.29 + 0.08 

12.92 + 0.94 

38.58/40 


The Cash statistic (Humphrey, Liu & Buote 2009) was used to fit the 
spectrum. 


of 1.72 mJy (Reinthal et al. 2012). Later, the flux level steadily de¬ 
creased. 

The long-term radio light curve is consistent with constant 
emission (probability of ~6.2xl0“^; ^^/d.o.f. = 61.25/25) with an 
average flux level of (0.136 ± 0.01) Jy. Considering only data taken 
during MAGIC observations (MJD = 55603 to MJD = 55622), the 
hypothesis of a non-variable emission at a mean flux level of 0.13 Jy 
has a probability of 75 per cent. Compared to 2010 November ob¬ 
servations, the radio data show a marginal flux increase from 2011 
mid-January to May, exceeding the mean flux level of the overall 
observation period. 

3.3 Variability study across wavebands 

To study the connection between variability patterns in individual 
wavebands, MWL data gathered during the MAGIC observations 


(from MJD 55568 to MJD 55622) were studied. Due to the sparse 
overlap of simultaneous multi-frequency data and the limited du¬ 
ration of the MWL observations, detailed methods such as e.g. 
calculating cross-correlation functions of simultaneous data sets 
were not applicable. We therefore settled for a simple linear 2D- 
regression analysis where we plotted the simultaneous data points 
from one waveband as a function of flux in the comparison wave¬ 
band. 

MWL data from five different wavebands (radio, optical, X- 
rays, HE y rays and VHE y rays) were available for comparison. 
The HE y-ray band was excluded due to the inability to detect lES 
0806+524 with Fermi-LAT on day timescales, as it was done in 
the other wavebands. The data in the optical. X-rays and VHE y 
rays were fairly simultaneous being on average less than 2h and 
individually no more than 3.25 h apart, while they all had an offset 
of roughly 12 h with respect to the radio observations. In the case 
of two radio measurements (MJD = 55616.4 and MJD = 55620.4) 
there were two observations in other wavebands, performed with a 
gap of ~12h before and after the corresponding radio observation. 
In those cases, we calculated the estimated flux level of the com¬ 
parison waveband at the time of the radio measurement assuming 
linear behavior of the flux between the two points. This led to the 
exclusion of radio to X-ray and radio to VHE data sets from the 
study due to a low number of concurrent data points. In total, we 
have four waveband pairs for the study; radio vs. optical, optical vs. 
X-rays, optical vs. VHE y rays and X-rays vs. VHE y rays. 

We settled on a 16 h window as a trade-off between strict si¬ 
multaneity and feasible daily overlap for selecting the MWL corre¬ 
lation points. In practice, this only affected the comparison of radio 
to other wavebands as all the rest of them had fairly simultaneous 
pointings. Nevertheless the number of available points was quite 
low ranging from three to eight with radio to X-rays and radio to 
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Figure 7. One-to-one connection plots derived for the MWL correlation studies cairied out between VHE y rays, X-rays as well as the optical and radio bands. 
The data are marked in black (filled circles), while the constant and the linear fits are represented by blue solid and red dashed lines respectively. 


Table 2. Results of the linear regression analysis. 



optical to 
radio 

optical to 
X-ray 

optical to 
VHE 

X-ray to 
VHE 

/ d.o.f (const, fit) 

5.8/5 

4.1/4 

19.3/7 

5.9/4 

X^ / d.o.f (lin. fit) 

3.0/4 

3.4/3 

17.8/6 

1.8/3 

fit likelihood (const, fit) 

0.33 

0.39 

0.0073 

0.21 

fit likelihood (lin. fit) 

0.56 

0.34 

0.0068 

0.62 


VHE having the lowest number of concurrent points. The result of 
this correlation study is shown in Fig 7. 

In order to determine whether there is a correlation between 
activities in different wavebands we compared the probability of a 
linear fit to the data to that of a constant fit (Table 2). 

We find that the linear fits are not significantly better than the 
constant fits. The case with the largest difference occurs for the X- 
ray to VHE flux, where a likelihood ratio test shows that the con¬ 
stant fit has a tail probability (computed using the Wilks theorem) 
of 0.04. 

It is worth noting that, while the VHE observations were trig¬ 
gered by the optical high state, there seems to be no one-to-one 
correlation between the fluxes in the two wavebands in the short 
term. Focussing on observations right before and after the flare, 
optical observations indicate a rather constant flux while the mea¬ 
surements performed in VHE y rays lead to the derivation of a 
short-term variability on a daily timescale. This suggests that the 
variability timescales are very different in these two bands and for 
future studies a longer time span should be used. 


4 SPECTRAL ENERGY DISTRIBUTION AND 
INTERPRETATION 

The SED of the source describing the high and low source state dur¬ 
ing the MAGIC observations, together with contemporaneous data 
from Fermi, Swift, the KVA and OVRO telescopes are presented in 
Fig. 8. EBL corrections have been applied to the VHE y-ray data 
using the model by Dominguez et al. (2011). 5wi/f/UVOT data and 
optical data in the R-band provided by the KVA telescope have 
been corrected for Galactic extinction and the host galaxy contribu¬ 
tion respectively according to Fitzpatrick (1999) and Nilsson et al. 
(2007). Simultaneous and quasi-simultaneous data have been com¬ 
bined accordingly to the high and low state in VHE y rays. In gen¬ 
eral, quasi-simultaneous data have to be considered when studying 
the broad-band variability of blazars because strictly simultaneous 
observations are not always available. In this particular study, by 
quasi-simultaneous observations we mean the usage of Swift (X- 
ray and UV) data from one day after the VHE flare detected by 
MAGIC. 

Unfortunately, the high state data do not include a simultane¬ 
ous HE y-ray detection, as Fermi did not significantly detect lES 
0806-1-524 during the flares in VHE y rays and X-rays. In order 
to provide simultaneous coverage in this energy band, the 95 per 
cent confidence level upper limit (300 MeV to 300 GeV) for 2011 
February 24 was included in the SED modeling (due to the poor 
statistics only one upper limit for the whole energy band has been 
calculated). In addition, an averaged SED from eight months of 
observations taken from 2010 November to 2011 June has been de¬ 
rived and is shown for comparison. Consequently, the latter data 
are not included in the modeling of the high state SED. For the low 
state of lES 0806-t-524, an averaged Fermi spectrum was produced 
considering data collected during the period from 2011 January to 
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March 2 and excluding the night of the VHE y-ray flare (February 
24) from the data sample. 

With respect to the Swift data, the high state comprises obser¬ 
vations carried out on February 25 since it is the closest X-ray spec¬ 
trum to the VHF flare detected by MAGIC on February 24, whereas 
the low state encompasses observations performed on February 1. 
Because of the high variability in the X-ray and VHF bands (es¬ 
pecially during the flaring state), the lack of simultaneity between 
these two energy bands places uncertainty in the interpretation of 
the broad band SFDs within the one-zone SSC model. Thus, we es¬ 
timate the X-ray spectrum for February 24 by scaling the spectrum 
from February 25 (Fig. 5) with a factor of 3.3, which is derived 
from the linear fit to the X-ray-to-TeV-flux (Fig. 7). Since we were 
not able to reproduce both the optical and VHF y-ray SFD points 
from previous MWL observations (Acciari et al. 2009), they are not 
shown. 

The MAGIC observations show a clear flux variability of 
about a factor of three in VHF y rays when comparing the high and 
low states (see Figs. 3 and 4). The Fermi-LAT light curve study in¬ 
dicated variability for the time period considered in this work (MJD 
= 55603 to MJD = 55622). 

The Swift/XRT spectrum from February 25 is marginally 
harder than that from February 1, and the X-ray flux above 2keV 
from February 25 is about three times higher than that from Febru¬ 
ary 1 (Table 1). The ratio between the assumed X-ray spectrum for 
February 24 (high state) and the flux of February f (low state) is 
about nine. 

The flux in the R-band started to increase well before the flar¬ 
ing activity in the VHF regime and in X-rays (see Fig. 5). It reached 
a maximum level at an almost constant flux with only marginal 
variations towards the end of February when the VHF y-ray flare 
occurred. Therefore, no clear variation between the high and low 
state SFDs is seen. While the inverse Compton peak indicates flux 
variability between both activity states, the synchrotron component 
shows only minor variation. Unfortunately, the weak detection in 
the HF y-ray band limits the determination of the inverse Compton 
peak. 

A one-zone SSC model is applied to reproduce the SFDs 
of both source states (for a detailed description of the model 
see Maraschi & Tavecchio 2003), where a spherical emission re¬ 
gion of radius R is assumed, filled with a tangled magnetic field of 
intensity B. A population of relativistic electrons is approximated 
by a smoothed broken power law that is parametrized by the min¬ 
imum y,ni„, break yt and maximum Lorentz factor as well as 
by the slopes iii and «2 before and after the break respectively. Rel¬ 
ativistic effects are taken into account by the Doppler factor d. The 
emission is self-absorbed at radio frequencies. This implies that the 
jet decelerates from the source region to the outer regions where the 
bulk of the radio emission originates. Therefore, radio data are not 
included in the SFD modelling. 

It is well known (Tavecchio, Maraschi, & Ghisellini 1998) 
that the one-zone SSC model is constrained once the basic SFD 
parameters and the variability timescale fyar - related to the source 
radius - are known. In the case studied here only an upper limit 
on the variability timescale is derived, fyar :£ 1 day. This timescale 
is much larger than the variability timescale that has been mea¬ 
sured during some flaring episodes in BL Lac objects, especially in 
the TeV band, for which sub-hour variability has been observed 
(Gaidos et al. 1996; Albert etal. 2007c; Aharonian et al. 2007). 
This uncertainty on Lar prevents us from strongly constraining the 
model parameters. In particular, the Doppler factor 6 and the mag¬ 
netic field B. We thus provide two possible realizations of the 


model, corresponding to two different values of the Doppler fac¬ 
tor (6 = 30 and <5 = 15) bracketing the typical range of values 
found in TeV BL Lac modeling (e.g. Tavecchio et al. 2010). 

The physical parameters derived reproducing the SFDs of 
both source states are reported in Table 3. To model the high state, 
we consider the quasi-simultaneous X-ray data from February 25 
observations [6 = 30) and the assumed simultaneous X-ray spec¬ 
trum (6 = 30 and 6 = 15). The values of the parameters are similar 
to those inferred for other HBL objects (see e.g. Tavecchio et al. 
2010 ). 

The differences between the high and the low state is mainly 
driven by the value of the Lorentz factor of the electrons dominat¬ 
ing the emission at the synchrotron peak, yj,, which almost doubles 
from the low to the high state, while the other parameters have al¬ 
most the same value (except for the electron normalization K that, 
for the case 6 - 30, halves from the low to the high state). As 
expected (Tavecchio, Maraschi, & Ghisellini 1998) the magnetic 
field intensity is inversely proportional to the Doppler factor, being 
therefore larger in the case d = 15. B field values of 0.05 - 0.1 G 
are commonly derived for HBLs (e.g. Tavecchio et al. 2010). 

We also report the power carried by the jet through electrons 
(Pc), magnetic field (Pb) and protons (Rp, derived assuming the 
presence of one cold proton per emitting electron). We note that 
while in FSRQ this choice is dictated by energetic requirements, 
in the case of BL Lac there are no stringent constraints on the 
composition (e.g. Celotti & Ghisellini 2008). However, this choice 
has a rather small impact on the derived total power, which is 
dominated by the leptonic component for both values of 6. The 
magnetic field and protons appear to contribute less to the jet 
power, which is mostly carried by the electrons. Its total value 
Pjct = Pe + Ps + Pf ~10‘*‘* erg s“* is also typical (e.g. Ghisellini et al. 
2010). While the jet power carried by the electrons and protons 
doubles during the source flare, the magnetic field strength is re¬ 
duced. 

The Doppler factors assumed for the SFD modelling re¬ 
ported in this paper (6 - 15 and 6 = 30) are larger 
than those typically derived from radio observations by about 
one order of magnitude. This discrepancy is commonly found 
in the TeV BL Lac population (Blasi et al. 2013). Such so- 
called Doppler crisis implies that the jet decelerates from the 
source region to the outer regions where the bulk of the ra¬ 
dio emission originates (e.g. Georganopoulos & Kazanas 2003; 
Ghisellini, Tavecchio & Chiaberge 2005). 

The assumed one-zone model, although applied here to the 
average SFD of the low state, is strictly valid only for snapshots, 
since the traveling and expanding blob is expected to change the 
emission properties because of adiabatic and radiative losses. As 
noted e.g. in Tagliaferri et al. (2008), the emission region may be a 
standing shock through which the jet plasma continuously passes. 

Assuming that the whole emission region Ren, is causally con¬ 
nected on the observed timescale, the physical parameters derive a 
minimal variability timescale of Lar.min = [Rem ■ (1 + z)](c ■ d) ~ 
0.3 (0.9) days for the high state assuming either the high (low) 
Doppler factor. The values found are both perfectly compatible 
with the variability timescale of one day inferred from the VHF 
light curve. 

The SSC model of the quasi-simultaneous data (Fig 8) im¬ 
plies an increase of the inverse Compton peak by a factor of three 
for the high state SFD with respect to the low state. When us¬ 
ing the assumed simultaneous X-ray spectrum, such increase is 
excluded. Considering the assumed simultaneous X-ray data, the 
SSC model indicates a shift of the synchrotron peak and the in- 
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Figure 8. SED of the high (left) and low (right) state of lES 0806+524 obtained during the MAGIC observations. The red (black) markers depict the 
simultaneous (quasi-simultaneous from February 25) MWL data. The MAGIC VHE y-ray data are shown with open squares, Fermi-LKi HE y-ray data with 
triangles and arrows for the upper limits, Swt/f-XRT data with filled squares, Swift-WCfl data with filled circles, KVA data with crosses and OVRO data 
with open diamonds. The solid (dotted) black curve depicts the one-zone SSC model matching the simultaneous (quasi-simultaneous) data assuming a high 
Doppler factor. We also show the modelling (dashed line) of the simultaneous data assuming a lower Doppler factor. For comparison purposes, the left panel 
shows also the fermf-LAT HE y-ray spectrum from 2010 November to 2011 lune (grey open triangles and grey arrows). See text for further details. 


Table 3. Input model parameters for the high and low state SEDs shown in Fig. 8. We report the minimum ymin, break yb and maximum Lorentz factor y,^iix 
and the low- and high-energy slope of the electron energy distribution, the magnetic field intensity, the electron density, the radius of the emitting region 
and its Doppler factor. We also report the derived power carried by electrons, magnetic field, and protons (assuming orre cold proton per emitting relativistic 
electron and the frequencies of the synchrotron and inverse Compton peak). (I): high state with quasi-simultaneous X-ray data; (Ila/b): high state with assumed 
simultaneous X-ray data modelled with a high/low Doppler factor; (Illa/b): low state modelled with a high/low Doppler factor. 


ymin 

[103] 

Tb 

[10^] 

Tmax 

[10®] 

rtl 

«2 

B 

[G] 

K 

[103 cm-3] 

R 

[10*® cm] 

6 

Re 

[10'*3 erg s-*] 

Pb 

[10^3 ergs-*] 

Pp 

[10'*3 ergs-*] 

"'^sync 

[Hz] 

VIC 

[Hz] 


1 

1.65 

7 

2 

3.85 

0.05 

19.0 

1.12 

30 

79.9 

0.1 

28.9 

10*5.26 

1023.84 

(I) 

1 

3.2 

7 

2 

3.85 

0.06 

3.0 

1.70 

30 

34.3 

0.35 

11 

10*5.97 

1024.®* 

(Ila) 

1 

3.6 

7 

2 

3.85 

0.10 

2.6 

3.00 

15 

24.0 

0.76 

7.3 

1015.90 

1024.48 

(Ilb) 

1 

2.1 

7 

2 

4.15 

0.05 

4.0 

1.70 

30 

43.0 

0.2 

14.3 

1015.33 

1023.87 

(Ilia) 

1 

2.8 

6 

2 

4.30 

0.10 

1.8 

3.20 

15 

18.6 

0.9 

5.8 

1015.56 

1024.09 

(Illb) 


verse Compton peak towards higher frequencies and a synchrotron 
dominance rather than an equal peak flux of both energy bumps. 
Such synchrotron dominance is more pronounced when using the 
high Doppler factor for the modelling. 

In comparison to the VHE data from Acciari et al. (2009), 
the inverse Compton peak is more constrained by the MAGIC 
data of the high source state. The comparison between the phys¬ 
ical parameters obtained from the low state SED observed by 
MAGIC to MWL observations of lES 0806+524 in 2008 March 
8 is not straightforward, as Acciari et al. (2009) applied an SSC 
model whose electron spectrum is approximated by an unbroken 
power law. The synchrotron and inverse Compton peak positions 
(vsync ~ 10*® Hz, Vic ~ lO^^'Hz, corresponding to ~41eV and 
~13 GeV) derived for the SED from 2008 (Acciari et al. 2009) and 
MAGIC observations presented in this paper are located at the same 
order of frequency. 

In contrast to Acciari et al. (2009), our SSC model implies a 
clear deviation from equipartition (with electrons dominating the 
total energy density), which is even more pronounced during the 
high state of the source. The origin of this difference with the 


VERITAS results is likely due to the different synchrotron peak 
frequency that, in the case of VERITAS, allowed a larger value of 
the magnetic field. 


5 SUMMARY & DISCUSSION 

In this article we report on the MAGIC observations of the HBL 
object lES 0806+524 from 2010 December to 2011 June, which 
were triggered by an optical high state. MAGIC detected a one-day 
long flare with a VHE flux three times larger than the flux during 
the low state. This flaring episode lasted less than one day. From 
this relatively short, high-activity state, whose observation in weak 
sources like lES 0806+524 is rather rare, a short-term variability of 
one-day timescale has been inferred. This in turn sets constraints on 
the size of the emitting region during the high VHE state. Excluding 
the flare night, the VHE y-ray flux from the source was in good 
agreement with the one observed by VERITAS in 2008, when the 
source was in a lower optical state. We present detailed spectra in 
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the energy range from a few hundred GeV to one TeV, describing 
both the high and the low VHE state. 

We studied the variability patterns of lES 0806+524 during 
the MAGIC VHE y-ray observations across the available wave¬ 
bands ranging from radio to VHE y rays in order to investigate the 
possible connection between the emission at different frequencies. 
The 2D linear regression we carried out did not reveal a connection 
between the flux levels in any of the possible waveband combina¬ 
tions. Although the VHE observations were optically triggered, no 
apparent evidence of a short-term correlation was found between 
these wavebands. 

As the multifrequency data cover both the high and the quies¬ 
cent VHE state activities of lES 0806+524, we performed one- 
zone SSC modelling of the SEDs obtained from the respective 
data sets assuming two different Doppler factors due to the un¬ 
certainty on the variability timescale. This model could adequately 
explain the broad-band emission during both source states using 
physical parameters similar to those from other HBL objects. Non¬ 
simultaneity of some of the MWL data hampered detailed studies 
on the broad-band variability of lES 0806+524. 
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